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We report the observation of harmonic generation and 
strong nonlinear coupling of two collective modes of a con- 
densed gas of rubidium atoms. Using a modified TOP trap we 
changed the trap anisotropy to a value where the frequency of 
the m — high-lying mode corresponds to twice the frequency 
of the m = low-lying mode, thus leading to strong nonlinear 
coupling between these modes. By changing the anisotropy 
of the trap and exciting the low-lying mode we observed sig- 
nificant frequency shifts of this fundamental mode and also 
the generation of its second harmonic. 

PACS numbers: 03.75.Fi, 05.30.Jp, 05.45.-a, 32.80.Pj 

Ever since the first Bose condensed gases were pro- 
duced 0, collective excitations have played a key role 
in the theoretical models and their experimental verifica- 
tion. Measurements of collective excitation frequencies 
may be compared with theoretical predictions and the 
first measurements of the lowest collective modes of a 
Bose gas were carried out soon after the first condensates 
were made [00|- Those experiments verified that the 
Gross-Pitaevskii (GP) equation, also called the nonlin- 
ear Schrodinger equation (NLSE), gives a very accurate 
prediction of the frequencies of the measured collective 
modes. The nonlinearity of the condensate is manifested 
as a term in the GP equation proportional to the mean- 
field or condensate density. In this experiment we have 
observed effects arising directly from this term. 

The spectroscopy of the excited states of the con- 
densate has been extended to include measurements of 
damping rates and frequency shifts as a function of tem- 
perature (in the range < T/T c < 1, where T c is the crit- 
ical temperature of condensation) B . The finite temper- 
ature measurements test dynamical aspects of the theory 
which require higher order terms in the calculations and 
as yet there is not full agreement || . 

Recent theoretical work has contributed to a better 
understanding of damping of the collective excitations at 
finite temperature |(| . Although mechanisms that lead to 
damping are complex one may consider two broad types 
of behavior. Firstly, energy may be lost from a mode 
by being transferred into a multitude of other modes 
through the nonlinear coupling that is intrinsic to con- 
densates. The second cause of damping is scattering of 
non-condensate particles in so-called Landau processes. 



The relative importance of these two processes depends 
on the shape of the condensate. In this work we have been 
able to vary the aspect ratio of the condensate around a 
value where nonlinear coupling is large and hence likely 
to dominate over other processes. 

In the hydrodynamic limit which applies to conden- 
sates with a large number of atoms, the GP equation 
reduces to the hydrodynamic equations of nondissipa- 
tive fluid dynamics. In this limit Stringari |t]] found 
analytic solutions for the frequencies of the collective 
modes and Edwards et al. J§| extended the calculations 
to the non-hydrodynamic region by direct integration of 
the GP equation. The measured oscillation frequencies 
correspond to those obtained from linear response the- 
ory only in the limit of small driving and response am- 
plitudes (zero amplitude limit). For strong driving the 
inherent nonlinearity in the condensate becomes appar- 
ent and the actual response frequency is shifted from the 
zero amplitude limit, loq. The shift Aco is predicted to be 
proportional to the square of the response amplitude to 
first order approximation, i.e. Alu/ujq = A 2 (5(A), where 
A is the excitation response amplitude, <5(A) is a nonlin- 
earity factor which varies for different trap geometries j^] 
characterized by A = oj z /oj r (the ratio of axial to radial 
trap frequencies) . Thus a frequency shift of the collective 
excitation can be achieved in two ways: either by increas- 
ing the driving and consequently the response amplitude 
of the oscillation or by changing the trap geometry and 
hence the nonlinearity factor. In previous work only the 
amplitude dependence of the frequency shift for a fixed 
trap geometry has been studied @-@]. This paper de- 
scribes the use of a trap with adjustable anisotropy JTq] 
to study the effects of trap geometry on collective exci- 
tations. 

The three lowest collective modes of a condensate in 
an axially symmetric trap (lj x = u> y — uj r ) are the m = 2 
mode, the low-lying m = mode and the high-lying 
7tl — mode, where m is the azimuthal angular momen- 
tum (the trivial centre of mass motion is not considered) . 
The m = 2 mode corresponds to a quadrupole type ex- 
citation in the radial plane. The low-lying m = mode 
corresponds to a radial oscillation of the width which 
is out of phase with an oscillation along the trap axis. 
The high-lying m = mode is an in-phase compressional 
mode along all directions (breathing mode). In our ex- 
periment we changed the trap anisotropy parameter A 
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around a resonance between the low-lying and the high- 
lying m — modes. At this resonance two quanta of 
excitation of the low-lying mode are converted into one 
quantum of excitation of the high-lying mode, i.e. the 
second harmonic of the low-lying collective mode is ex- 
cited. The theoretical plot (for the hydrodynamic limit) 
in Fig.|l]a shows that at a trap anisotropy of A w 1.95 
(resonance) the frequency of the high-lying m = mode 
corresponds to exactly twice the frequency of the low- 
lying mode. 

The nonlinearity factor 8(X) also changes dramatically 
across the resonance. A theoretical plot of S(X) against 
A is given in Fig.|l]b using the formula given in 

Our experimental apparatus for creating Bose conden- 
sates is described elsewhere 11 1 and we only briefly men- 
tion the relevant features here. Atoms are loaded into a 
TOP trap and then evaporatively cooled to the quantum 
degenerate regime. We cool the atoms to well below the 
critical temperature where no thermal cloud component 
is observable (~ 0.5T C ). For this experiment the usual 
TOP trap configuration with a rotating field in the radial 
plane was modified by the addition of a pair of 'z-coils' 
in Helmholtz configuration which produced an oscillating 
axial bias field along the (vertical) quadrupole coil axis. 
The total field in the trap is 



B(t) = B' q (xe x + ye y - 2ze z ) + 
B r (cos LOtt e x + sinwti e y ) + B z cos2ct> t t e 2 



(1) 



The first term is the static quadrupole field and is written 
in terms of its gradient B' in the radial direction. The 
second term is the conventional TOP field of magnitude 
B r rotating in the xy-plane at a frequency u>t ■ The third 
term is the additional z-bias field of magnitude B z mod- 
ulated at twice the frequency of the field in the xy-plane. 
Hence the locus of the quadrupole field does not simply 
follow a planar circular path like that of a standard TOP 
trap but it also moves up and down following a saddle 
shape. 

The time average of this field configuration yields the 
value for a trap frequency ratio that is smaller than 
A = \/8 in a TOP trap with B z = 0. By increasing 
B z the anisotropy A is tuned continuously from 2.83 to 
1.6. The practical limit on A in our present apparatus 
arises because of the noise created by the amplifier driv- 
ing the z-coils which induces Zeeman substate changing 
transitions for large amplitudes of B z . This limit may 
be overcome to create a spherical trap (A = 1) as de- 
scribed in , which also gives more details of the tech- 
nique. We calibrated the trap by measuring the centre 
of mass oscillations of a thermal cloud for various ampli- 
tudes of B z . The measured frequencies agree to better 
than 1% with a theoretical prediction obtained by nu- 
merically time-averaging the field in Eq.[l], i.e. integrat- 
ing over one complete TOP cycle. Note that the radial 
trap frequency remains approximately constant whereas 



the axial frequency decreases as the amplitude B z in- 
creases. The trap frequencies for this experiment were 
126 Hz radially and the axial frequencies varied from 356 
Hz at B z /B r = (A = V8) to 194 Hz at B z /B r = 2.3 
(A = 1.6). 

Once the trap had been calibrated we investigated the 
behavior of the m — mode for various anisotropics. To 
excite the m = low-lying mode, the TOP-field ampli- 
tude B r was modulated sinusoidally at a frequency of 225 
Hz (matching the hydrodynamic mode frequency in a 126 
Hz radial trap with B z = 0) for a period of 8 cycles (35 
ms). Note that the actual measured mode frequency for 
the above trap parameters differs from 225 Hz because 
the finite number of atoms leads to deviations from the 
hydrodynamic limit. Also, the change of mode frequency 
with trap anisotropy (see Fig.|l|a) means that the driv- 
ing frequency does not exactly match the hydrodynamic 
mode frequency for traps with anisotropics different from 
2.83 (the value for B z = 0), but it is close enough to ex- 
cite the fundamental mode for all our measurements. 

After being excited the condensate was left to oscillate 
for variable times of between 2 to 30 ms before the mag- 
netic trap was switched off and the condensate allowed to 
expand freely for 12 ms to make a time-of-flight (TOF) 
measurement. Then absorption images were taken of the 
condensate from which we extracted both the radial and 
axial widths and the total number of atoms in the conden- 
sate (typically 2 x 10 4 ). The radial TOP field produced 
an 8% modulation of the trap spring constant and we 
measured response amplitudes of around 20% in TOF. 
These correspond to even smaller amplitudes of oscilla- 
tion for the condensate before expansion which produce 
no measurable frequency shift in the ordinary TOP trap 
(for A = 2.83 the nonlinear coupling 5 for the m = mode 
is close to zero, see Fig. 0b). Even for strong driving and 
large response amplitudes no significant frequency shift 
of the m = mode from the zero amplitude limit has 
been observed before 0. 

The anisotropy of our trap was tuned across the res- 
onance between the low and the high-lying m = 
mode. We observed the occurrence of the second har- 
monic frequency when we tuned our trap close to res- 
onance (Fig.||b) and on resonance it was the dominant 
component (Fig.||c). It arises from the excitation due to 
nonlinear coupling of the m = high-lying mode and 
could only be observed in the oscillations of the axial 
width. The radial width oscillated at a single frequency 
corresponding to the m = low-lying mode and did not 
have a second frequency component (Fig.^d,c). This is 
also what we found in simulations based on the hydrody- 
namic equations and indicates that the geometries of the 
low and high-lying modes are such that the excitation 
of the high-lying mode can only be observed by measur- 
ing the axial widths. We obtained the low- lying mode 
frequency from a single frequency fit to the radial os- 
cillations and also from a two frequency fit to the axial 
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oscillations and found good agreement between the two 
values (exponential decay factors were included in both 
fits). The frequency of the second harmonic (correspond- 
ing to the high-lying mode) was obtained as the second 
frequency component in the fit along the axial direction. 
We found that the fundamental frequency was strongly 
suppressed on resonance in both the radial and axial os- 
cillations as shown in Fig.||c,f. This indicates that the 
low-lying mode is not populated anymore and all the en- 
ergy of the excitation has been transferred to the high- 
lying mode ]l3fl . The large error for the point at reso- 
nance (A = 1.93) in Fig.|| reflects the difficulty in obtain- 
ing the fundamental component from the data in Fig.||c. 

The frequencies of the fundamental mode for various 
trap geometries were normalized by the corresponding 
radial trap frequency and plotted as a function of the 
trap anisotropy parameter A in Fig.||a. 

The dotted line is the prediction of the small amplitude 
mode frequencies in the hydrodynamic limit as given in 
§. The solid line was calculated by Hutchinson juj from 
the GP equation for the atomic number and trap frequen- 
cies of our experiment. Figure |^a shows that the mea- 
sured mode frequency is above the hydrodynamic predic- 
tion for the ordinary TOP trap with A = 2.83 in agree- 
ment with previous experiments 0]. Far from resonance 
the data agree very well with the finite number predic- 
tion (solid line). However, when approaching the reso- 
nance from above the nonlinear coupling becomes large 
enough to shift the mode frequency below the predicted 
curve. On the lower (A < 1.95) side of the resonance 
the frequency shift becomes positive in agreement with 
the predicition illustrated in Fig.|l|b. One can see that the 
measured frequency shifts follow a dispersive curve which 
crosses zero at A = 1.93 ± 0.02, in good agreement with 
the prediction of A = 1.95 for the hydrodynamic regime 
Pj . Figure |^b shows a plot of the ratio of the squared 
amplitudes (higher to lower-lying mode) which is pro- 
portional to the ratio of the mode energies. Away from 
resonance we observed no second harmonic contribution 
and the ratio (and its error) is taken as zero. These 
data have been fitted by a simple Lorentzian centered at 
A = 1.94 ± 0.02 with HWHM= 0.035 ±0.005. The damp- 
ing rates are plotted against the anisotropy A in Fig. ^c. 
They were around 20s -1 for the fundamental and 50s -1 
for the second harmonic frequency. 

In conclusion, we have a TOP trap in which the 
anisotropy can be tuned whilst keeping the trap frequen- 
cies high |l(J. This offers the opportunity to investigate 
a range of interesting phenomena in Bose-Einstein con- 
densates related to the trap geometry, in particular the 
greatly enhanced nonlinear effects for certain anisotropy 
parameters A. We observed the generation of the sec- 
ond harmonics of the m = low-lying mode frequency as 
well as a dispersive resonance in the mode frequency at 
A = 1.93 ±0.02. Future improvements and modifications 
of our experiment will allow us to study another predicted 



resonance at A = 1.5 as well as additional nonlinear phe- 
nomena such as the collapse and revival of oscillations 
and the onset of chaos for stronger driving amplitude 
fll5| . A particularly interesting case is the spherical trap 
(A = 1) where all frequencies are degenerate. Quantita- 
tive theoretical predictions have been made for the damp- 
ing rates in a spherical trap Jl6[ which can be compared 
with the experiment to provide a stringent test. 
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FIG. 1. (a) The coupling between the low and the 
high-lying mode is strongest for a trap anisotropy of A = 1.95. 
In this geometry the frequency of the high-lying mode is twice 
the frequency of the low-lying mode, (b) The nonlinearity 
factor 5(A) is plotted against the anisotropy A. In this region 
strong nonlinear features are observed around the resonance 
for second harmonic generation. Both plots are based on the 
hydrodynamic theory. 
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FIG. 2. Shape oscillations of the BEC are plotted against 
time for various geometries. Sub-plots (a,b,c) show oscilla- 
tions in the axial width and sub-plots (d,e,f) oscillations in 
the corresponding radial width, (a) Away from resonance 
(A = 2.64) only a single frequency appears, (b) Close to reso- 
nance (A = 1.97) there is a strong admixture of a second har- 
monic contribution, (c) On resonance (A = 1.93) the second 
harmonic contribution dominates. (d,e) The radial widths 
oscillate at the frequency of the fundamental mode (m = 
low-lying mode) for A = 2.64, 1.97. (f) The oscillation of the 
fundamental mode disappears on resonance. 
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FIG. 3. (a) The frequency of the m = low-lying mode in 
units of the radial trap frequency as a function of the trap 
anisotropy for various geometries. The dotted line is the pre- 
diction of the hydrodynamic theory and the solid line was 
calculated from the GP equation, (b) The measured energy 
ratio of the higher to the lower-lying mode is plotted against 
the trap anisotropy. The solid line is a Lorentzian fit to the 
data, (c) The damping rates for the fundamental (black cir- 
cles) and the second harmonic (white circles) frequency are 
plotted against the anisotropy A. 
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